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! ABSTRACT 

o 

^SJ ■ Although ACDM model is very successful in many aspects, it has been seriously challenged. 

> ^ ' Recently, warm dark matter (WDM) remarkably rose as an alternative of cold dark matter (CDM). 

I In the literature, many attempts have been made to determine the equation-of-state parameter 

(EoS) of WDM. However, in most of the previous works, it is usually assumed that the EoS of 
dark matter (DM) is constant (and usually the EoS of dark energy is also constant). Obviously, 
this assumption is fairly restrictive. It is more natural to assume a variable EoS for WDM (and 
dark energy). In the present work, we try to constrain the EoS of variable WDM with the current 
cosmological observations. We find that the best fits indicate WDM, while CDM is still consistent 
with the current observational data. However, ACDM is still better than WDM models from the 
C_J , viewpoint of goodness-of-fit. So, in order to distinguish WDM and CDM, the further observations 

■ on the small/galactic scale are required. On the other hand, in this work we also consider WDM 

• ' whose EoS is constant, while the role of dark energy is played by various models. We find that the 

i-C ' cosmological constraint on the constant EoS of WDM is fairly robust. 
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I. INTRODUCTION 

Nowadays, dark energy has become one of the most active fields in physics and astronomy, since the 
great discovery of the current accelerated expansion of our universe in 1998 [J, the other hand, 

dark matter (DM) was invoked to interpret the rotation curves of spiral galaxies for many years 
Since it is well known that hot dark matter (HDM) cannot be competent for the cosmological structure 
formation, cold dark matter (CDM) has become the leading candidate. In fact, the well-known ACDM 
model has been established as the standard model in cosmology today [ij . 

Although ACDM model is very successful in many aspects, recently it has been seriously challenged. 
According to the brief reviews in e.g. 0), these serious challenges to ACDM model include, for example, 
(1) ACDM predicts significantly smaller amplitude and scale of large-scale velocity flows than observa- 
tions; (2) ACDM predicts fainter Type la supernova (SNIa) at high redshift z; (3) ACDM predicts more 
dwarf or irregular galaxies in voids than observed; (4) ACDM predicts shallow low concentration and 
density profiles of cluster haloes in contrast to observations; (5) ACDM predicts galaxy halo mass profiles 
with cuspy cores and low outer density while observations indicate a central core of constant density and 
a flatfish high dark mass density outer proflle; (6) ACDM predicts a smaller fraction of disk galaxies 
due to recent mergers expected to disrupt cold rotationally supported disks. Even when one replaces the 
cosmological constant A with other (dynamical) dark energy candidates, these challenges still cannot be 
successfully addressed. In particular, the main source of the challenges on the small/galactic scale might 
be CDM. We refer to e.g. [3| for details. 

On the other hand, recently warm dark matter (WDM) remarkably rose as an alternative of CDM. 
The leading WDM candidates are the keV scale sterile neutrinos. In fact, the keV scale WDM is an 
intermediate case between the eV scale HDM and the GeV scale CDM. Unlike CDM which is challenged 
on the small/galactic scale (as mentioned above), it is claimed that WDM can successfully reproduce 
the astronomical observations over all the scales (from small/galactic to large/cosmological scales) 
The key is the connection between the mass of DM particles and the free-streaming length (structure 
smaller than £[s will be erased). The eV scale HDM is too light and hence all structures below Mpc scale 
will be erased; the GeV scale CDM is too heavy and hence the structures below kpc scale cannot be 
erased (therefore CDM is challenged on the small/galactic scale); in between, the keV scale WDM works 
well We refer to e.g. [3] for a comprehensive review. 

It is well known that the equation-of-state parameter (EoS) plays an important role in cosmology. In 
particular, the EoS of CDM and radiation/HDM are and 1/3, respectively. In between, the EoS of 
WDM, Wm, satisfles < Wm ^1/3 obviously. Of course, the realistic value of Wm should be determined 
from the astronomical observations. A non-zero Wm indicates WDM, rather than CDM. In the literature, 
many attempts have been made to determine the EoS of DM. For example, in [H, assuming a constant 
Wm and the role of dark energy is played by a cosmological constant A whose EoS is —1 exactly, it is found 
that -1.50 X 10"^ < u;™ < 1.13 X 10"^ if DM produces no entropy, and -8.78 x 10"^ < < 1.86 x 10"^ 
if the adiabatic sound speed vanishes. Note that in ^ the allowed range of the EoS of DM was relaxed and 
could be negative, but this is somewhat unnatural and hence it is not the case which will be investigated 
in the present work. Following the method proposed in Q which suggested to measure the EoS of 
DM by combining kinematic and gravitational lensing data, the authors of found that the EoS of 
DM is consistent with pressureless dark matter within the errors. The authors of Q considered the 
observational constraints on a cosmological model with variable EoS of matter and dark energy, namely, 
Wm = l/[3(a;" + 1)] and Wde = wx"/ (a;" -I- 1), where x = o/a* with a* being a reference value of the scale 
factor a, and a > is a constant model parameter. The model considered in Q was proposed to unify 
the radiation dominated phase and the dark energy dominated phase. The matter behaves as radiation 
when a; 1 and DM when a; ^ 1. Note that this motivation is not for WDM, and its EoS Wm, Wde 
are ad hoc in some sense. The authors of @ considered the cosmological constraints on WDM whose 
EoS is a constant, while the EoS of dark energy is also a constant. They claimed that the cosmological 
data favor Wm = 0.006 ± 0.001 (suggesting WDM), and Wde = —1.11 ± 0.03 (corresponding to phantom 
dark energy). The CDM whose EoS Wm — and the cosmological constant A whose EoS Wde = — 1 are 
disfavored beyond 3(t confidence level. 

Note that in most of the previous works, it is usually assumed that the EoS of DM is constant (and 
usually the EoS of dark energy is also constant). Obviously, this assumption is fairly restrictive. It is 
more natural to assume a variable EoS for WDM (and dark energy). In fact, this is the main subject 
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of the present work. We will try to constrain the EoS of variable WDM with the current cosmological 
observations. In Sec. |IT1 we briefly introduce the observational data used in this work. In Sec. IIII[ we 
consider the cosmological constraints on WDM whose EoS is variable. In fact, we adopt the familiar 
Chevallier-Polarski-Linder (CPL) parameterization for WDM, namely, Wm — Wmo + Wma{^ — a). Unlike 
the somewhat ad hoc parameterization in noting that the Taylor series expansion of any function 
F{x) is given by F{x) = F{xo) + Fi {x - xq) + {F2/ 2!) {x - Xof + (F3/ 3!) {x - xof + . . . , the CPL 
parameterization for WDM can be regarded as the Taylor series expansion of Wm with respect to the 
scale factor a up to first order (linear expansion), and hence it is naturally motivated. In this section, we 
consider three cases, namely, the role of dark energy is played by a cosmological constant A, and dark 
energy described by constant EoS, CPL parameterized EoS, respectively. In Sec. IIVI noting that in most 
of the previous works constant EoS of both WDM and dark energy are assumed, here we also consider 
WDM whose EoS is constant, but dark energy is a cosmological constant A, or dark energy whose EoS 
are constant and CPL parameterized, respectively. We try to see whether the cosmological constraint on 
the constant EoS of WDM is robust for various types of dark energy, especially when the EoS of dark 
energy is variable. In Sec. |Vl some concluding remarks are given. 



II. OBSERVATIONAL DATA 

Recently, the Supernova Cosmology Project (SCP) Collaboration released the Union2.1 compilation 
which consists of 580 Type la supernovae (SNIa) [l^. The Union2.1 compilation is the largest published 
an d sp ectroscopically confirmed SNIa sample to date. The data points of the 580 Union2.1 SNIa compiled 
in are given in terms of the distance modulus iXobs{zi)- On the other hand, the theoretical distance 
modulus is defined as 

lJ.th{zi) = ^logiQ D L[zi) + ^jio , (1) 
where fiQ = 42.38 — blogi^ h and h is the Hubble constant Hq in units of 100 km/s/Mpc, whereas 

in which p denotes the model parameters; z is the redshift; E = H/Hq, in which H = a/a is the 
Hubble parameter; a dot denotes the derivative with respect to cosmic time t; a = (1 + z)^^ is the scale 
factor (we have set oq = 1; the subscript "0" indicates the present value of corresponding quantity). 
Correspondingly, the from 580 Union2.1 SNIa is given by 

where a is the corresponding la error. The parameter /iq is a nuisance parameter but it is independent 
of the data points. One can perform a uniform marginalization over fj,Q. However, there is an alternative 
way. Following [ill . IT^ , the minimization with respect to hq can be made by expanding the xfi of Eq. Q 
with respect to /j,o as 



X^(p) = i-2/ioS + MgC, (4) 



where 



Af ^ s:^[fJ.obs{zi)- iJ,th{zi;iJ.o = 0,p)] 



Of \ sr^ f^obsiz.^) - fith{zf,^io = 0,p) 
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Eq. (|4]) has a minimum for ^iq — B/C aX 

x» = i(p)-:^. (5) 

Since x^. min = X^i. min (up to a constant) obviously, we can instead minimize which is independent of 
jUo- In addition to SNIa, the other useful observations include the cosmic microwave background (CMB) 
anisotropy [l^ and the large-scale structure (LSS) However, using the full data of CMB and LSS 

to perform a global fitting consumes a large amount of computation time and power. As an alternative, 
one can instead use the shift parameter R from CMB, and the distance parameter A of the measurement 
of the baryon acoustic oscillation (BAO) peak in the distribution of SDSS luminous red galaxies. In the 
literature, the shift parameter R and the distance parameter A have been used extensively. It is argued in 
e.g. [IHl that they are model-independent and contain the main information of the observations of CMB 
and BAO, respectively. As is well known, the shift parameter R of CMB is defined by [isl [l6j 

where the redshift of recombination = 1091.3 which was determined by the Wilkinson Microwave 
Anisotropy Probe 7- year (WMAP7) data and Q.mO = ^t^G pmo / {^H^) is the present fractional density 
of matter. The shift parameter R relates the angular diameter distance to the last scattering surface, the 
comoving size of the sound horizon at and the angular scale of the first acoustic peak in CMB power 
spectrum of temperature fluctuations [l^, [l^. The value of R has been determined to be 1.725 ± 0.018 
from the WMAP7 data [l3|- On the other hand, the distance parameter A of the measurement of the 
BAO peak in the distribution of SDSS luminous red galaxies 14] is given by 

A ^ ^,ilE{z,)-^/^ 

where zi, = 0.35. In [l3|, the value of A has been determined to be 0.469 (ns/0.98)"°-3^ ± 0.017. Here the 
scalar spectral index Ug is taken to be 0.963, which comes from the WMAP7 data [13]. So, the total x^ 
is given by 

2 ~2 I 2 I 2 /o\ 

X = + XCMB + Xbao ) (8) 

where x^ is given in Eq. ©, Xcmb = - Robsf/(T\ and x\ao = " Aobsfl(T\- The best-fit 
model parameters are determined by minimizing the total x^- As in (l8l - [20l |. the 68.3% confidence level 
is determined by Ax^ = X^ - xlrtn < 1-0, 2.3, 3.53, 4.72 and 5.89 for Up = 1, 2, 3, 4 and 5 respectively, 
where Up is the number of free model parameters. Similarly, the 95.4% confidence level is determined by 
Ax^ = X^ - Xm»« < 4-0, 6.18, 8.02, 9.72 and 11.31 for Up = 1, 2, 3, 4 and 5, respectively. 
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III. COSMOLOGICAL CONSTRAINTS ON VARIABLE WDM 



In this section, we consider the cosmological constraints on WDM whose EoS is variable. Throughout 
this work, we consider a flat Friedmann-Robertson- Walker (FRW) universe which contains dark energy 
and WDM (we assume that radiation and the baryon component can be ignored). Here, we adopt the 
familiar CPL parameterization for WDM, namely, its EoS is given by 

Wm = Wma + Wmai^ " a) , (9) 

where Wmo and Wma are both constants. Unlike the somewhat ad hoc parameterization used in as 
mentioned in Sec. HI the CPL parameterization for WDM can be regarded as the Taylor series expansion of 
Wm with respect to the scale factor a up to first order (linear expansion), and hence it is naturally moti- 
vated. To ensure < Wm < 1/3 in the whole history (0 < a < 1) as mentioned in Sec. HI we require that 

0<Wmo<l/3, < w„io + w,„Q < 1/3 . (10) 

In this section, we consider three cases, namely, the role of dark energy is played by a cosmological 
constant A, and dark energy described by constant EoS, CPL parameterized EoS, respectively. 
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FIG. 1: The 68.3% and 95.4% confidence level contours in the Wmo — Wma plane for the AVWDM model. The 
best-fit parameters are also indicated by a black solid point. 




FIG. 2: The 68.3% and 95.4% confidence level contours in the WmO — Wma plane for the XVWDM model. The 
best-fit parameters are also indicated by a black solid point. 
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At first, we consider the case in which the role of dark energy is played by a cosmological constant A. 
We call it AVWDM model. As is well known, the corresponding E{z) is given by [l^ [21I, 



E{z) 



exp 



l + z 



-,1/2 



(1 - n. 



(11) 



There are 3 free parameters in this model. By minimizing the corresponding total in Eq- ®i we find 
the best-fit parameters = 0.2774, Wmo = 0.0 and Wma = 0.0036, while x^m = 562.227. In Fig. [U 
we present the corresponding 68.3% and 95.4% confidence level contours in the Wmo — Wma plane for the 
AVWDM model. It is easy to see that the best fit indicates WDM, while CDM is still consistent with 
the current observational data. Note that for the best fit, Wm = at z = 0, namely it is CDM today 
while it is WDM in the past (a < 1). 

Next, we consider the case in which the EoS of dark energy Wx is a constant. We call it XVWDM 
model. In this case, the corresponding E{z) is given by fiol [2l[ [22| 



E{z) 



exp 



3u>„ 



l + z 



+ (1 - f7„,o) (1 + 



1/2 



(12) 



There are 4 free parameters in this model. By minimizing the corresponding total x^ in Eq. (0), we 
find the best-fit parameters fimo = 0.2773, Wmo = 0.0072, Wma = —0.0071, and = —1.0072, while 
Xmin — 562.225. In Fig. [51 we present the corresponding 68.3% and 95.4% confidence level contours in 
the Wmo — Wma plane for the XVWDM model. Again, it is easy to see that the best fit indicates WDM, 
while CDM is still consistent with the current observational data. Note that for the best fit, it is always 
WDM in the whole history (0 < a < 1), and Wm — Wmo — 0.0072 today. Comparing Figs. [T]and[51 we 
can easily find that the confidence level contours become larger. 

Finally, we consider the case in which the EoS of dark energy Wy^ is also variable. Here, we adopt CPL 
parameterization for dark energy, namely 



= WxO + «^xa(l - a) , 



(13) 
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where Wxo and w^„ , are both constants. We call it CVWDM model. In this case, the corresponding E{z) 



There are 5 free parameters in this model. By minimizing the corresponding total in Eq. ([5]), we find 
the best-fit parameters Qmo — 0.2773, w^o — 0.0072, w,na = —0.0072, w^o = —1.0081, and Wxa = 0.0062, 
while Xmin — 562.225. In Fig.[3l we present the corresponding 68.3% and 95.4% confidence level contours 
in the Wmo — Wma plane and the Wxo — Wxa plane for the CVWDM model. From Fig. |31 we see that ACDM 
is still consistent with the current observational data. Note that for the best fit, it is always WDM in the 
whole history (0 < a < 1), and Wm ~ WmO ~ 0.0072 today. 



FIG. 4: The 68.3% and 95.4% confidence level contours in the Q,„m - Wm plane for the AWDM model. The 
best-fit parameters are also indicated by a black solid point. 



IV. COSMOLOGICAL CONSTRAINTS ON WDM WITH CONSTANT EOS 

Noting that in most of the previous works constant EoS of both WDM and dark energy are assumed, 
here we also consider WDM whose EoS Wm is constant, but dark energy is a cosmological constant A, 
or dark energy whose EoS are constant and CPL parameterized, respectively. We try to see whether the 
cosmological constraint on the constant EoS of WDM is robust for various types of dark energy, especially 
when the EoS of dark energy is variable. 

We firstly consider the case with a cosmological constant A. We call it AWDM model. Note that the 
EoS of WDM should satisfy < Wm < 1/3 as mentioned in Sec. HI As is well known, the corresponding 



is given by [T^, 





(14) 



0.014 





E{z)= r!™o(l + ^)'^'+"''"' + (l-f^mo) 



nl/2 



(15) 
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There are 2 free parameters in this modeL By minimizing the corresponding total in Eq. ([8]), we find 
the best-fit parameters Slmo — 0.2770 and Wm = 0.0023, while Xmin = 562.228. In Fig. [H we present 
the corresponding 68.3% and 95.4% confidence level contours in the fimo — w,n plane for the AWDM 
model. It is easy to see that the best fit indicates WDM, while CDM is still consistent with the current 
observational data. 

Then, we turn to the case in which the EoS of dark energy Wx is a constant. We call it XWDM model. 
In this case, the corresponding E{z) is given by [lol . [2ll [2^ 

1/2 



E{z) = r!„o(i + + (1 - n^o) (1 + 2)3(1+-") 



(16) 



There are 3 free parameters in this model. By minimizing the corresponding total in Eq. ([8]), we 
find the best-fit parameters flmo = 0.2776, Wm = 0.0025, and Wx = —1.0035, while Xmin = 562.225. In 
Fig. [SJ we present the corresponding 68.3% and 95.4% confidence level contours in the Wm — w^ plane for 
the XWDM model. Again, the best fit indicates WDM, while CDM is still consistent with the current 
observational data. 
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FIG. 5: The 68.3% and 95.4% confidence level contours in the w„ 
parameters are also indicated by a black solid point. 



- Wx plane for the XWDM model. The best-fit 



Finally, we consider the case with CPL parameterized dark energy, namely, the EoS of dark e nergy is 
given by Eq. (fT3)) . We call it CWDM model. In this case, the corresponding E{z) is given by |19l l2lL |23| 

1/2 



a„o(l + z)3(i+-'") + (1 - r!,„o) (1 + 



exp 



3wx 



1 



(17) 



There are 4 free parameters in this model. By minimizing the corresponding total in Eq. ([S]), we 
find the best-fit parameters f^^o = 0.2773, Wm = 0.0023, w^o = -1.0074, and w^a = 0.0299, while 
Xmin = 562.224. In Fig. |6l we present the corresponding 68.3% and 95.4% confidence level contours in 
the Wm — Wxo plane and the Wm — Wxa plane for the CWDM model. From Fig. [51 we see that ACDM is 
still consistent with the current observational data. 

Comparing these three cases, it is easy to find that their best-fit Wm are almost the same value. This 
indicates that the cosmological constraint on the constant EoS of WDM is fairly robust, namely it is 
insensitive to the dark energy models. 
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FIG. 6: The 68.3% and 95.4% confidence level contours in the Wm — WxO plane and the Wm — Wxm plane for the 
CWDM model. The best-fit parameters are also indicated by a black solid point. 

V. CONCLUDING REMARKS 

Although ACDM model is very successful in many aspects, it has been seriously challenged. Recently, 
warm dark matter (WDM) remarkably rose as an alternative of cold dark matter (CDM). In the literature, 
many attempts have been made to determine the equation-of-state parameter (EoS) of WDM. However, 
in most of the previous works, it is usually assumed that the EoS of dark matter (DM) is constant (and 
usually the EoS of dark energy is also constant). Obviously, this assumption is fairly restrictive. It is 
more natural to assume a variable EoS for WDM (and dark energy). In the present work, we try to 
constrain the EoS of variable WDM with the current cosmological observations. We find that the best 
fits indicate WDM, while CDM is still consistent with the current observational data. On the other hand, 
in this work we also consider WDM whose EoS is constant, while the role of dark energy is played by 
various models. We find that the cosmological constraint on the constant EoS of WDM is fairly robust. 

As mentioned above, in the six cases considered in this work, all the best fits indicate WDM, while 
CDM is still consistent with the current observational data. So, it is worthwhile to compare these WDM 
models with ACDM. To this end, we also fit ACDM model to the same observational data, and find 
the best-fit fimo = 0.2736, while Xmin = 562.546. Following 0, here we adopt three criterions 
used extensively in the literature, namely Xmin/'^^f^ Bayesian Information Criterion (BIC) and Akaike 
Information Criterion (AIC). Note that the degree of freedom dof ^ N — k, whereas N and k are the 
number of data points and the number of free model parameters, respectively. The BIC is defined by [23l | 

BIC = -2\nCrr,ax+klllN , (18) 

where C„iax is the maximum likelihood. In the Gaussian cases, x^j„ = —2h\ Cmax- Thus, the difference 
in BIC between any two models is given by 

ABIC = AxL„ + Afcln7V. (19) 

On the other hand, the AIC is defined by [11] 

KlC^-2\nCraax+2k. (20) 

Correspondingly, the difference in AIC between any two models reads 



AAIC = AxLn + 2Afc. 



(21) 
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Model ACDM AVWDM XVWDM CVWDM AWDM XWDM CWDM 

xLin 562.546 562.227 562.225 562.225 562.228 562.225 562.224 

k 1345234 

■)dmr,ldof 0.968238 0.971031 0.972708 0.974393 0.969359 0.971028 0.972706 

ABIC 12.4139 18.7784 25.1449 6.04847 12.4119 18.7774 

AAIC 3.681 5.679 7.679 1.682 3.679 5.678 

Rank 1 4 6 7 2 3 5 



TABLE I: Comparing all the six WDM models with ACDM. 



In Table H we present Xmm/^o/> ABIC and AAIC for ACDM and all the six WDM models considered 
in this work. Notice that ACDM has been chosen to be the fiducial model when we calculate ABIC and 
AAIC. From Table HI we see that the rank of models is coincident in all the three criterions (Xmi„/rfo/, 
BIC and AIC). Obviously, ACDM model is the best one. In summary, although the best fits to the 
cosmological observations (SNIa, CMB and BAO) indicate WDM, we cannot say WDM is favored, since 
ACDM is still better than WDM models from the viewpoint of Xmi„/do/, BIC and AIC. So, in order to 
distinguish WDM and CDM, the further observations on the small/galactic scale are required. 
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